Effective theory for low-energy nuclear energy density functionals 
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We introduce a new class of effective interactions to be used within the energy-density-functional 
approaches. They are based on regularized zero-range interactions and constitute a consistent 
application of the effective-theory methodology to low-energy phenomena in nuclei. They allow 
for defining the order of expansion in terms of the order of derivatives acting on the finite-range 
potential. Numerical calculations show a rapid convergence of the expansion and independence of 
results of the regularization scale. 
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A consistent formulation of the low-energy energy- 
density-functional (EDF) approach in terms of effective 
theory is long overdue in nuclear physics. Effective the- 
ories are now used in different domains of physics with a 
tremendous success. They quantify and formalize basic 
physical intuitions about aspects of the description that 
are essential, and those that can be treated in terms of 
controlled and correctible approximations. One spectac- 
ular example relates to derivations of nuclear forces based 
on the chiral perturbation theory augmented by the 
effective-field-theory treatment of corrective terms [ll, Q ■ 
Analogous effective approach can also be applied in a 
much simpler context of quantal one-body dynamics 0] . 

Successful derivations from chiral effective field the- 
ory of relativistic [1, [H and nonrelativistic 0] nu- 
clear EDFs were recently accomplished. They lead to 
local EDFs with density-dependent couplings associated 
with the underlying pion-exchange interaction. These 
approaches strive to map out the in-medium nucleonic 
effects at the two-pion-mass (2TO7r ~ 260MeV/c^) or 
Fermi-momentum {kp ~ 260MeV/ftc ~ 1.3 fm~^) scales. 
However, the energy scale of low-energy nuclear phenom- 
ena in finite nuclei is much lower. For example, to dis- 
sociate a nucleon from a nucleus one has to increase its 
kinetic energy by (JE'kin — fi^kpSk/M ~ 0.25hc6k, that 
is, a typical separation energy of (5-Ekin — 8 MeV corre- 
sponds to the momentum increase of 6k ~ 32MeV /he. 
Moreover, the low-energy EDFs are supposed to describe 
nuclear excitations and shell-effects at even lower ener- 
gies of 1 MeV or 6k ~ AMeV/hc and below. 

Therefore, what in the QCD-driven chiral dynamics 
is the small-energy scale, becomes a short-range high- 
energy scale of nucleon-nucleon force acting on relatively 
weakly bound nucleons in nuclei. The question of bridg- 
ing these two different energy scales is the subject of 
current intense studies in nuclear physics Q. In finite 
nuclei, surface effects decrease the infinite-matter parti- 
cle binding energies by about a factor of two. In the 
present study, we aim at building a phenomenological 
low-energy EDFs, with gradient terms providing for the 



surface terms, corresponding the nuclear leptodermous 
expansion (ol. [lOj. 

Ideas of the effective theory have recently been applied 
in constructing the EDFs based on hig her-order deriva- 
tives of nuclear one-body densities [11|- When analyzed 
in terms of the density-matrix expansion [l^ . the ob- 
tained quasilocal EDFs converge fast with the order of 
derivatives [l^, which is a very encouraging feature of 
the approach. The quasilocal EDFs are also obtained 
by the Hartree-Fock averaging of higher-order zero-range 
pseudopotentials [T^ . 

However, one aspect of the approach developed so far is 
still missing, namely, the quasilocal EDFs and zero-range 
pseudopotentials do not provide us with any expansion 
scale, against which the higher-order derivatives could be 
classified. In the chiral effective field theory, such a scale 
is given by the cut-off energy in the momentum-space 
integrals, and allows for formulating a consistent power- 
counting scheme of the diagrammatic expansion |15| . In 
low-energy nuclear physics, an analogous cut-off regu- 
larization of the zero-range interactions has already been 
proposed Il6l, | l7l , and the dimensional regularization was 
also studied [l8|, however, a consistent expansion scheme 
is not yet available. 

In the present study, we introduce the expansion scale 
by employing the regularized higher-order zero-range 
(pseudo)potentials - exactly as it was illustrated in the 
simple examples discussed in Ref. [H. Among a plethora 
of various ways to regularize delta interactions, we con- 
sider the one based on the Gaussian function, 



S(r) = lim Qair) = lim 5-. 
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Then, the resulting central two-body regularized pseu- 
dopotential reads, 

F(r;,r^;ri,r2)= (2) 
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where k = ^^(Vi - V2) and k' = Yii^'i ^ ^2) ^-rc the 
standard relative-momentum operators, and the Wigner, 
Bartlett, Heisenberg, and Majorana terms are given by 
the standard spin and isospin exchange operators, Pi = 
1, P2 = Pa, P3 = -Pr, Pa = -PaPr- The relativc- 
momentum operators are the building blocks of the scalar 
higher-order derivative operators Oi{k,k'). Generaliza- 
tion of the pseudopotcntial to spin-orbit and tensor terms 
are straightforward [l4| , and in the present study are not 
discussed. 

In our definition of the pseudopotential ©, operators 
Oi{k,k') act on the product of the locality delta func- 
tions S{r[ — ri)S{r'2 — and regularized delta function 
ga{ri — ■r2). Derivatives of the delta functions have to be 
understood in the usual sense of distributions, namely, 
when the pseudopotential ^ is inserted into a matrix 
element, or when it is averaged with respect to den- 
sities, the integration by parts transfers all derivatives 
onto wave functions or densities. After the integration 
by parts, the locality deltas reduce the integral to two 
variables, ri and r2, only. 

To give a specific example, up to the second-order, that 
is, up to the next-to-leading-order (NLO) expansion, op- 
erators Oi{k,k') read 



1, 



(i) 

where are the channel-dependent coupling constants. 
In the limit of a — > 0, Eq. ([3]) links the regularized pseu- 
dopotential ([2]) to the standard zero-range Skyrme inter- 
action (l9| . 

Before considering the general pseudopotcntial ([2]) , let 
us first assume that the differential operators Oi{k,k') 
depend only on the sum of relative momenta, Oi [k -\- k'). 
At NLO, this corresponds to equal second-order coupling 
constants, = . Such particular differential opera- 
tors commute with the locality deltas 5{r'i — ri)5{r2—r2), 
and thus can be applied directly onto the regularized 
delta Qaifi — r2). In such a case, pseudopotential ^ 
reduces to a simple local potential 



k'^)+Tl'h-k\ 



(3) 
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where Vi{r) = Oi{k)ga{r) are functions of the relative 
distance r = ri — r2- Moreover, since Oi{k) are scalar 
differential operators, the potentials must have forms of 
power series of Laplacians A in r, that is. 



(5) 



where Vj^,*' are the coupling constants at order 2n. 

Within the spirit of the effective theory, the series in 
Eq. ([5]) is cut at a given rimax- Then, at a given fixed scale 
a, coupling constant ¥2^ have to be adjusted to data. If 



the expansion is relevant, the results should only weakly 
depend on the regularization scale a, provided it is nei- 
ther too large, when the expansion scale intrudes into the 
physical scale of low-energy phenomena, nor too small, 
where unresolved high-energy features are incorrectly de- 
scribed. At any given order one can pick a specific scale 
a where the results agree best with the data. 

Before such a complete program is carried out, in the 
present exploratory study we use the fact that a success- 
ful parametrization of the local potential already exists in 
the form of the Gogny interaction [2l[ , whose finite- 
range part can be written as 



1=1 1=1 A:=l,2 



(6) 



where Gaussians with two ranges, ai = 0.7 fm and 
02 = 1.2 fm, have been employed. Therefore, instead 
of adjusting the regularized potential to data, we first 
attempt to derive it directly from the Gogny interaction. 
In the present study, we focus on the finite-range part of 
the Gogny interaction and we keep its zero-range spin- 
orbit and density-dependent terms unchanged. We then 
use the results obtained for the original Gogny interac- 
tion as metadata, wherewith we compare those obtained 
for the derived regularized potentials. 

At NLO (riniax — 1), the eight coupling constants, Vq^^ 
and V2^\ of the local regularized potential @ can be 
easily obtained from the eight parameters of the central 
Gogny force, G\ and Gj . Indeed, at this order, poten- 
tials ([5|) are combinations of a Gaussian and Gaussian 
multiplied by r^, which span similar sets of functions as 
those of two Gaussians with two different ranges. 

Below we use a prescription that can be applied at 
any order nmax, and which is based on global character- 
istics of the potentials. Namely, we determine coupling 
constants VJ^n in Eq. ([5]) by requiring that the lowest 
moments of both potentials are equal, that is. 



-'^^2m — 
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for m — 0,1, ... , rimax- This conditions gives the coupling 
constants of the regularized potential in simple analytical 
forms. 
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fc=l,2 



where G^*'' and are the parameters of the Gogny in- 
teraction, as in Eq. ([5]). 

To determine properties of finite nuclei, in the present 
study we used the EDF approach for regularized poten- 
tial defined in Eq. , along with the coupling constants 
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^ derived for the DIS Gogny force [2l|. Numerical cal- 
culations were performed by implementing regularized 
potentials in the code HFODD (v2.53r) [21 [11, which 
was adapted to solving self-consistent equations for finite- 
range interactions. 

In Fig. [Tl we present results of calculations performed 
for eight doubly magic nuclei: ^^O, ^^Ca, "^^Ca, ^^Ni, 
■^^Ni, i°°Sn, ^^"^Sn, and 2°®Pb. In each case, we show 
relative deviations of binding energies (left panels) and 
rms radii (right panels) as compared to values obtained 
with the original Gogny interaction. One can see that 
with increasing order of expansion, the convergence is 
very rapid; indeed, for every next order, the obtained 
relative deviations decrease by about a factor of four, 
and at N'^LO they are well below 1%. 
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for other standard functionals used within the EDF ap- 
proaches. Nevertheless, even without such extensive ad- 
justments, basic features of the obtained effective theory 
can already be illustrated by studying scaling properties 
of the results. To this end, in Fig. [2] we show deviations 
of binding energies and radii relative to those obtained 
for ^°^Pb. One can see that the results scaled in this way 
are almost independent of the regularization scale, which 
is an extremely encouraging result, consistent with the 
effective-theory principles. 
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FIG. 2: (Color online) Same as in Fig. [1] but for the devi- 
ations relative to those obtained for ^"^Pb. For clarity, near 
a ~ 1 fm, where the '^''^Pb normalization values cross zero 
(cf. Fig. [T]), several points were removed from the plot. 



FIG. 1: (Color online) Deviations (in percent) of the EDF 
binding energies [left panels (a), (c), and (e)] and rms radii 
[right panels (b), (d), and (f)], calculated for the regularized 
potential relative to values obtained with the original 
Gogny interaction. Top, middle, and bottom panels corre- 
spond to the NLO, N^LO, and N^LO expansions, that is, to 
2nniax=2, 4, and 6, respectively. The legend gives the atomic 
masses of doubly magic nuclei listed in the text. At NLO, 
missing points correspond to large values of a, where con- 
verged self-consistent solutions could not be obtained. 

We note here that the results presented in Fig. [T] were 
obtained without any adjustment of parameters. How- 
ever, one can expect that such small adjustments can 
easily correct for the smooth trends observed in Fig. [1] 
This is so because binding energies of nuclei, which result 
form cancellations of the kinetic and potential energies, 
are extremely sensitive to parameters of the interaction. 
For example, we checked that modifications of the Gogny 
density-dependent term of the order of 0.05% are able 
to bring the N^LO mass of ^°*Pb exactly to the Gogny 
value. 

The regularized potentials proposed in this study (with 
or without the restricted dependence on the relative mo- 
menta introduced above), will later be used in extensive 
optimization procedures - similarly as it has been done 



Effective theories are built to provide us with expan- 
sions that converge better when the bound-state energy 
or scattering energy is lower. Spectacular examples of 
such a convergence pattern exist for weakly bound elec- 
tron states in quantum mechanics of a one-body Coulomb 
problem Q or for phase shifts of the two-body nuclcon- 
nucleon scattering . Convergence of an effective theory 
can be best visualized by the so-called Lepage plot 0], 
which shows the dependence of the error in the descrip- 
tion of a given observable on energy. Such error should 
decrease with decreasing energy and with increasing or- 
der of expansion. 

In many physical systems, scale of energy is inverted 
with respect to the scale of distance, that is, small dis- 
tances correspond to large energies and vice versa. In 
finite nuclei, the situation is more complex, because the 
average energy per particle and average internucleon dis- 
tance are almost constant throughout the mass chart, 
which is a simple consequence of the saturation of nu- 
clear forces and near-constancy of density inside nuclei. 
Should then the effective theories work better in light 
(small) nuclei or in heavy (large) ones? 

Good arguments in favor of the second option are pro- 
vided by the density- matrix expansion (DME) [T2j . which 
uses the fact that the range of nuclear forces is smaller 
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than characteristic distances at which densities in nuclei 
vary, see Refs. [H, [l^ and references cited therein, and 
by the ideas of the leptodermous expansion [1, ■ In- 
deed, since the largest variations of densities occur at 
nuclear surfaces, effective theories should work best in 
large nuclei, where the surface region is relatively small 
with respect to the nuclear bulk. 

To study this problem, in Fig. [3] we show reduced de- 
viations of binding energies, Ab — AB/A, and radii, 
Ar = AR/A-^/^, in function of the particle number A. 
This allows us to compare results obtained for differ- 
ent nuclei in the same scale. Since the relative devia- 
tions are almost independent of the regularization scale 
a (cf. Fig. [2]), here we plot them at a = 0.85 fm. 
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FIG. 3: (Color online) Same as in Fig.[2l but for the deviations 
of binding energies scaled by the numbers of particles, Ab = 
AB/A, and those of radii scaled by the nuclear sizes, Ar — 
AR/A^^^. Note that the N^LO results almost perfectly hide 
behind those obtained at N^LO. 

We see that the results obtained for nuclei beyond 
about A = 4:8 scale with mass as those for ^'^^Pb. From 
the point of view of the effective theory and proposed 
expansion, one can interpret this result as a high simi- 
larity of these nuclei, that is, their sizes are not different 
enough to induce significant differences in their conver- 
gence properties, and they all would appear at the same 
point of the Lepage plot. We also see that properties of 
lighter nuclei, with A < 48, converge slightly better than 
those of the heavier ones. However, only a future com- 
plete analysis, performed for the entire mass tabic and 
for optimized functionals, can here provide conclusive re- 
sults. 

We complete our analysis by discussing values of cou- 
pling constants ([H]) in the so-called natural units, see 
Refs. [2^ [2^ and references cited therein. Since the 
EDF generated by our regularized potential contains only 
terms quadratic in densities, in the natural units the cou- 
pling constants have the following dimensionlcss values. 
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where is an overall normalization factor and is a 
characteristic scale of the derivative operator A. 

Coupling constants (0) are plotted in Fig.|31 Their 
almost-linear n-dependence, which we see in the logarith- 
mic scale, allows us to derive the slope of this dependence 
that corresponds to the value of A ~ TOOMeV/Sc ~ 
3.5 fm~^ (dashed line). For our regularized poten- 
tials at N'^LO, we thus obtain a similar scale A as 



that characteri zing the zero-range interactions at NLO 
(2nniax=2) [2^ l26l | Next, we fix the normalization fac- 
tor / so as to make the values of dimensionlcss cou- 
pling constants Wj'^ of the order of unity, which gives 
/ ~ 35MeV/(fic)3/2 ~ 1/(77 MeVi/2fm3/2). As we see, 
this value is significantly smaller than the pion decay con- 
stant /tt ~ 93McY/{hc)^/^, which normalizes the EDFs 
derived directly from the chiral effective field theory [2^ . 
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-S FIG. 4: (Color online) Coupling constants l/jn derived at 



0.85 fm from the Gogny interaction, as dictated by Eq. ((Sjl. 
The dashed line illustrates the slope of the n-dependence of 
the coupling constants, plotted here for A ~ 700MeV/fic. 



Finally, for the derived values of A and /, in Fig. [5] 
we show the dependence of the dimensionlcss coupling 
constants on the regularization scale a. Wc sec that 
for n ~ 1, 2, and 3 they have values of the order of unity 
at all scales, although the fourth-order coupling constants 
{n = 2) cross zero around a = 0.95 fm. The zero-order 
coupling constants (n = 0) appear to be less natural, with 
values of the order of 0.1 (cf. also Fig. |4|). However, we 
should keep in mind that in the present study we derived 
the coupling constants from the Gogny interaction, which 
has a predefined scale of the order of 1 fm. Therefore, we 
can expect that the future adjustments of the coupling 
constants directly to data may lead to even weaker scale- 
dependence. 
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FIG. 5: (Color online) Coupling constants in natural units 
as functions of the regularization scale a, plotted for n — 
0-3 in panels (a)-(d). 

In summary, in this study we introduced a new class of 
energy density functionals that are based on pseudopo- 
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tentials related to regularized delta interactions. We used 
the regularization scheme employing the standard Gaus- 
sian form, which introduces a regularization scale a cor- 
responding to the range of the interaction. Different or- 
ders of expansion then correspond to different orders of 
derivatives used in the pseudopotential. 

In our opinion, future prospects for using the proposed 
regularized (pseudo)potentials are high. First, similarly 
as in the simple one-body examples [3| , they may present 
better convergence properties than similar expansions 
based on the zero-range interactions. Second, they allow 
for convergent summations of contributions from high 
single-particle momenta, which is not the case for zero- 
range interactions. And third, they allow for formulating 



a consistent expansion in terms of the orders of deriva- 
tives, with the convergence properties gauged against the 
regularization scale. This last feature is unique among all 
the EDF approaches based on zero-range and finite-range 
interactions developed so far, and gives us a potential of 
building an order-by-order correctible theory. 
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